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Abstract 

The utility of prompt photon plus associated jet production for constraining 
the size and shape of the polarized gluon density of the proton AG is examined 
at = 40 GeV, appropriate for the proposed HERA-N polarized collider 
experiment. The calculation is performed at next-to-leading order {0{aa1)) 
in QCD. The reliability of the predictions are examined in some detail. 
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I. 



Polarized deep inelastic scattering (DIS) experiments have proved invaluable for pro- 
viding information on the spin- dependent parton densities of the nucleon through 
measurement of the polarized structure functions such as gi{x, Q"^). Data from such experi- 
ments have traditionally been used in phenomenological analyses which attempt to constrain 
the polarized parton densities and evolve them in through the use of the spin- dependent 
Altarelli-Parisi evolution equations. Up until fairly recently, all such analyses had to be 
performed in leading order (LO) QCD since the spin-dependent two loop splitting functions 
were unknown. Since the two- loop splitting functions have been calculated @J^, recent anal- 
yses have been performed in next-to-leading order (NLO) QCD and parametrizations 
of the polarized parton densities have been provided which can be used in calculations of 
other polarized processes. 

In these analyses extensive use is made of theoretical constraints such as the Bjorken 
sum rule in order to extract the parton densities from the polarized structure function 
measurements. Despite this the data does not provide tight constraints on sizes, and even 
more so the shapes of either the polarized quark (Ag) or gluon (AG) densities. This is true 
in particular for AG. Each group provides more than one different parametrization of the 
parton densities, each set fitting the structure function data, but they can have very different 
individual parton densities reflecting the limited constraints on these quantities provided by 
DIS experiments. It is thus very important, if our understanding of the spin structure of the 
nucleon is to be improved, that accurate measurements of these distributions be obtained 
from other processes. 

Since the cross section for prompt photon production is dominated by the subprocess 
qg 7g in hadronic collisions already in leading order, it has proved very useful for pro- 
viding information on the unpolarized gluon densities, g{x,Q'^), of hadrons at fixed target 
facilities. It has therefore been suggested that it may prove equally useful in pinning down 
the polarized gluon densities In this context it has been examined quite extensively. 
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The NLO corrections were calculated in and [|TT|, where numerical estimates were also 



presented and it was established that the LO results were perturbatively stable. In both 
those analyses only LO polarized parton densities were available and the fragmentation con- 
tributions were not included. More recently, prompt photon production with polarized beam 
and target was examined in ref. [T^JT^ using NLO polarized parton densities and including 
the fragmentation contributions. 

Although single inclusive prompt photon production will definitely be very important for 
constraining the size of AG, information on the detailed x-shape of the distribution may not 
be as easily extracted. This is because the calculation of the inclusive cross section involves 
one convolution over the momentum fractions, x, of the initial partons. The practical effect 
of this is that a measurement of the kinematic variables of the photon is not sufficient to 
determine the value of x. If, on the other hand, one or more of the jets produced in the 
reaction is also tagged, no convolution is involved in the calculation and the cross section is 
directly proportional to the parton densities. Thus the double longitudinal asymmetry All 
is directly proportional to the ratio AG/g in kinematic regions where other subprocesses 
such as qq scattering can be neglected. 

Schematically the cross section is given in terms of the polarized and unpolarized hard 
subprocess cross sections (A)o"ij and parton densities (A)/f(xa,M^) (the presence of the A 
before a quantity means the polarized version of the quantity) by 



(A)a = Y.i^)fnxa, M') * (A)/|(x,, M') * 



The * indicates either a product or convolution and explicit dependence of the hard sub- 
process cross sections on the kinematic variables of the observed products and the renor- 
malization and factorization scales have been omitted. The sum over i,j extends over all 
initial partons and DJ{z,Mp) is the unpolarized fragmentation function for a parton, c, 
into a photon at scale Mj^. Thus, just as in the single photon case, there are two classes of 
contributions to the cross section labelled the direct and fragmentation contributions. 
The polarized and unpolarized parton densities are defined by 



AMx, M') = f+{x, M') - f-{x, M') 



(1.2) 



and 

Mx, M') = f+(x, M') + fr(x, M'), (1.3) 

where the and /~ represent the distribution of partons of type i with positive and negative 
hehcities respectively, with respect to that of the parent hadron. The hard subprocess cross 
sections are defined by 

Aa,, = I - (1.4) 

and 

a,, = i (<7++ + <7+-) . (1.5) 

One of the main quantities studied in polarized experiments is the longitudinal asymmetry 
All. This quantity gives a measure of the sensitivity of the process to polarization effects. 
The double spin asymmetry studied in this paper is defined by 

All = — , (1.6) 
a 

the ratio of the cross section for longitudinal polarization of the incoming hadrons to the 
corresponding unpolarized cross section. 

In LO the direct processes contributing to the cross section are qg — > 7g and qq — > ^yg. 
In a LO jet calculation, the simple approximation 'parton = jet' is used irrespective of the 
jet definition used in the experiments, hence either the final state gluon or quark will form 
the jet. It is only in NLO that more than one parton may make up a jet, and one finds a 
dependence of the cross section on the jet definition and parameters. In addition there are 
the fragmentation processes: 

qg qg 
qq qq 
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qq' qq' 
qq qq 

QQ qg 

qq -> gg 
99^ 99 

99 m (1-7) 

where one of the final state partons fragments to produce the photon, i.e., q{g) — > 7 + X. 
The cross section of interest here is the triple differential cross section 



dpTpdrf dr]^ ' 

where 77'^ and t]^ are the pseudorapidities of the photon and jet respectively and pip is the 
transverse momentum of the photon. The direct contribution to the cross section is given 
in LO by 

2nplj:^aAfnx,,M')x,Af^(xj,M')^. (1.8) 



dpldT^^dr^J -'^1 ^-a^^^y-^^-- J-o^^3y-3^- > ^£ 

The corresponding fragmentation contributions will involve a convolution over the fragmen- 
tation variable z. s — XaX^S where is the cms energy in the hadron-hadron system and 
as usual, t = (pi — P'yY where pi and are the momenta of one of the initial state partons 
and the observed final state trigger photon respectively. The Bjorken variables, Xa,&, are 
given in terms of the kinematic variables of the photon and jet by 



7 
Pt 



^ ^e-^-" + e-^') . (1.9) 

In leading order pj. — p^ and thus a measurement of pj^, rf and 77"^ means that Xa,h can 
be determined. In other words, a measurement of d^Aa'^^ /dpj^drfdr]^ corresponds to a 
measurement of d^Aa"''^/dxadxb- In NLO due to the presence of a third parton in the 
process this simple formula is no longer applicable. In that case, in order to determine Xa,b 



unambiguously, one would need to measure the kinematic variables of all three final state 
partons. The values of Xa^b are only approximately determined by the transverse momenta 
and pseudorapidities of the photon and jet in this case, due to the need to implement a jet 
definition. If photon isolation requirements are also implemented, then the determination 
of Xa^b from the photon and jet kinematic variables would be even less precise. Such a 
measurement would still nevertheless yield some very useful information on the x-dependence 
of the polarized parton densities as will be shown. 

In NLO, O(aaf), there are virtual corrections to the LO non-fragmentation processes, 
as well as the further three-body processes: 



Q +9 - 


q + 9 + 1 


(1.10a) 


9+9- 


^ q + q + 'y 


(1.10b) 


q +q- 


^9+9+1 


(1.10c) 


q +q- 


-> g + g + 7 


(l.lOd) 


q +q- 


g + g + 7 


(l.lOe) 


q + q- 


^ g' + g' + 7 


(l.lOf) 


q + q' - 


-> g' + g + 7 


(l.lOg) 



The virtual corrections as well as the three-body matrix elements were calculated in ref. |TT 
and I use these matrix elements in this calculation. 

In principle the the fragmentation processes of eq.(1.7) should now be calculated to 0{al) 
and convoluted with the NLO photon fragmentation functions whose leading behavior is 
0{a/as), but the hard subprocess matrix elements are not yet available in the polarized case. 
Hence unless otherwise stated, in both the polarized and unpolarized cases, I include the 
leading order contributions to these processes only. Numerically the fragmentation processes 
are not significant except at extremely low pt due to the low cms energy of the fixed target 
experiment as I shall show later, but for a theoretically consistent calculation they should 
nevertheless be included as they help to reduce uncertainties from scale dependence. 
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For this calculation I use the Monte Carlo method first used in for the unpolarized 
case of the photon plus jet calculation performed here, as well as for calculations of some 
other important processes, to deal with the phase space integrals. In this letter I shall not 
give any details of the calculation except to say that the calculation was performed in the 



MS scheme using the t'Hooft-Veltman (HVBM) scheme |[T6| to treat 75. I refer to ref. |T7 
for the details of the calculation. At NLO, since there are two-to-three scattering processes 
included in the calculation it is possible that two partons may be too close together to be 
resolved as two separate jets. In this case a jet definition is required. In this calculation I 
use the cone definition proposed at Snowmass [|l^, which defines a jet as hadronic energy 



deposited in a cone of radius Rj = y (Arjy + (A^)^. At the partonic level, if one parton 
forms the jet then the kinematic variable of the parton is set equal to that of the jet, pj = Pi, 
7] J = rji and (pj = (pi, and the jet is centered on the parton. If another parton falls inside the 
jet cone then the jet variables are the weighted averages of those of the two partons: 

Pj = Pi+ P2 

Vj = — {PiVi + P2V2) 
Pj 

(pJ = — {Pl(pl + P2(p2) ■ (1-11) 

Pj 

All results are displayed for pp collisions at the center-of-mass energy y/s = 40 GeV 



appropriate for the proposed HERA-N experiment at DESY |Tj] . For the unpolarized cross 
section the CTEQ3M parton densities are used throughout, and the value of corre- 
sponding to this distribution is also used. Use of other unpolarized parton densities at the 
x-values probed here do not yield significantly different results. For the polarized case the 
GRSV and GS distributions are used with the corresponding values for A-jjg. The 
authors of ref. and |^ have proposed various parametrizations of the polarized parton 
densities differing mainly in the choice of input for the polarized gluon density AG. In the 
case of the GRSV distributions we use the 'valence' set which corresponds to a fit of the 
available DIS data (referred to by the authors as the 'fitted' AG scenario), the large gluon 
fit which assumes that AG(x, Qq) = g{x, Qq) at input (the 'AG = (7' scenario) and the small 



gluon fit which uses AG{x,Qq) = at the input scale (the 'AG = 0' scenario), which in 
this case starts at the very low value of Ql = 0.34 GeV^. The latter two distributions are 
intended to represent extreme choices for AG. These parametrizations give gluon densities 
which differ in their absolute sizes as well as in their x-shape. The GS parametrizations 
provide three fits to the data; GS A, GS B and GS C. It has been shown that the GS A and 
GS B distributions do not differ very much from the AG = g and fitted AG sets of GRSV 
respectively, whereas the the GS C set is widely different from any of the others. I shall 
present distributions using the three GRSV sets discussed above, along with the GS C set 
for comparison. For the fragmentation functions I use the LO asymptotic parametrization of 
ref. |jl9|. As will be shown, the choice of fragmentation functions makes very little difference 
to the predictions, since these processes account for only a small fraction of the cross section. 

The renormalization, factorization, and fragmentation scales are set to a common value 
fi = Pt unless otherwise stated. Dependence on fi is examined in one of the figures below. 
Since there are two 'particles' in the final state, the jet and the photon, both of whose 
transverse momenta are large, an alternative choice might be = pj. or some function of 
Pj- and p^. The results of the calculations show that the magnitudes of pj. and p^ tend 
to be comparable and that dependence of the asymmetries on /i is slight, although the 
individual cross sections may vary significantly with /i. Therefore, choices of fi different 
from fi = Pj' should not produce significantly different predictions for the asymmetries . 
The two loop expression for as{Q'^) is used throughout, with the number of fiavors fixed at 
Nf = 4, although the contribution from charm was verified to be negligible at the energies 
considered. The value of the jet cone size is fixed at Rj = 0.5 unless otherwise stated. 

Fig. la shows the triple differential cross section as a function of pj^ of the photon for the 
various parametrizations. The unpolarized cross section is shown for comparison. The curves 
were obtained by averaging over Arj = 1 and Ap^ = 1 GeV and the restriction pj. > pj. was 
imposed. The upper dotted curve is for the GS A parametrization, verifying that it is very 
similar in both shape and size to the AG = g parametrization of GRSV. All the remaining 
parametrizations give distributions which are distinctly different in both their shapes and 
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sizes. If one compares these distributions with the corresponding curves for single inclusive 
prompt photon production |12| one finds, apart from the expected fall in the absolute size of 
the cross section, that the corresponding distributions differ also in their shapes. This is most 
obvious for the GS C parametrization which has the most distinctly different shape. Fig. lb 
also shows significantly larger asymmetries than the inclusive photon case for equivalent 
p^'s. Most of these differences can be traced to the fact that for the inclusive photon case, 
a given value of pj. corresponds to a much less sharply defined value of x than in the photon 
plus jet case. In Fig.la,b the range 4 < < 10 GeV corresponds to 0.2 < Xa,b < 0.5. It 
was also verified that for all parametrizations, the qg scattering process accounts for more 
than 90% of the cross sections. Also plotted in fig. lb are the projected statistical errors in 
the asymmetries, 6 All, as estimated in ref. from the formula 



SAll = 0.17 /^a{pb). (1.12) 
In WM SAll has been calculated by assuming an integrated luminosity of 240 pb^^ and 



beam and target polarizations of Pb = 0.6 and Pt = 0.8, and including a trigger and 
reconstruction efficiency of 50% with no acceptance correction. Similar errors for 6 All were 
estimated in |T2[ by integrating the unpolarized cross section over bins of sizes Arj = 1 and 



ApT = 1 GeV. We use the same proceedure here. The results indicate that for < 7 — 8 
GeV the asymmetries will be distinguishable. 

Fig.lc and Fig. Id are similar to Figs.la,b but in this case both the jet and photon are 
restricted to the forward rapidity regions, t]'^ and rj'^ are centered at +1, with At] = 0.2. At 

= 5 GeV this corresponds to probing Xa,b in bins centered at the points Xa,b = 0.07 and 
Xb,a = 0.55. The asymmetries are smaller than in fig. lb but again there are clear distinctions 
in sizes and shapes for the different parametrizations. 

In Fig. 2 we look at distributions in r]'^ at = 5 GeV and allow p'^ to vary between 5 
and 20 GeV. In Fig. 2a the jet is restricted to be in the central region, -0.5 < < 0.5, 
and in Fig. 2c it is restricted to the forward region, 0.5 <ri^< 1.5. The asymmetry plots of 
Fig. 2b shows important differences between the various parametrizations but that of Fig. 2d 

9 



is very striking. The asymmetries are rather flat, although showing differences in shape for 
the various parametrizations in the positive rapidity regions, but are very large and increase 
very sharply in the negative rapidity region, almost approaching 1 at the edge of phase space 
for the AG = g scenario. 

An examination of Fig. 2c shows that this effect originates with the difference in shapes 
between the polarized and the unpolarized rapidity distributions. It was argued in ref. ||2^ 



that positive rapidity correlations at collider energies are an inherent property of the hard 
subprocess matrix elements. Thus if one restricts the jet to be at positive rapidity then the 
rapidity distribution of the photon would be expected to peak in the positive rapidity region. 
This was verifled recently in ref. |2l[] where prompt photon plus charm quark correlations 



were studied at NLO. Clearly in Fig. 2c this expectation is conflrmed for the unpolarized 
cross section, but in contrast the polarized cross sections show negative correlations. Since 
the propagators in the hard subprocess matrix elements are the same this can only be due to 
the differences between the polarized and unpolarized parton densities. This fact was verifled 
by using artiflcial polarized parton densities of exactly the same shape as the unpolarized 
ones along with the polarized matrix elements. In that case positive rapidity correlations 
between the photon and the jet were obtained. 

In ref. |2^ it is suggested that it is the behavior at small-Xa, Xh of the product 
XaXbf'^{xa,Q'^)f''{xb,Q^) along with the structure of the hard subprocess matrix elements 
which generate the positive rapidity correlations between two flnal state particles. It is well 
established that the small-x behavior of the ratio AG/g is ~ x as x — > 0. Hence there is an 
additional power of x in the polarized parton distributions at small-x. This also applies in 
the set for which AG = g at the input scale. At the scale we are considering {Q"^ ~ 25 GeV^) 
the evolution equations determine the shape of the parton densities regardless of the shape 
at the input scale. Therefore given the fact that the behavior of the ratio AG/g at small-x 
is generated by the polarized Altarelli-Parisi equations themselves and that the GRSV 
input scale is so low {Qq = 0.34 GeV^), it is understandable that even the AG = g scenario 
will also show the same small-x behavior characteristic of the polarized parton densities. It 
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is easy to show from eq.(1.9) pO[ that when only two final state particles are present, as is 
the case in LO 

XaXb = ^^il + cosh{r]^ -ri^)). (1.13) 

The effect of this additional factor is to suppress the polarized rapidity distributions at the 
point 7]'^ = 7]^ and produce two symmetrical peaks on both sides of this point. There are 
many other factors such as, for example, available phase phase which also affect the rapidity 
correlations and may tend to obscure the effect of the small-x behavior, but from Fig. 2d 
the effect is a significant rise in the predicted asymmetries at negative rapidities. This 
strong sensitivity to polarization effects in this region should serve as a very good test of the 
underlying QCD mechanism as well as a check on assumptions about the small-x behavior 
of the polarized parton densities. 

In order to test the sensitivity of the predictions to the choice of factorization and renor- 
malization scales. Fig. 3a shows predictions for the asymmetry of a sample of the the pj. 
distributions of Fig. la for three different scales. All hard scales are kept equal and varied 
simultaneously. Although, as in the case of single prompt photon production, varying the 
scales over the full range shown can change the individual cross sections by up to 50%, the 
predicted asymmetries are relatively stable. In Fig. 3b the sensitivity of the predicted asym- 
metries to the cone size of the jet Rj and to the inclusion of the fragmentation contributions 
is tested. The solid lines are again the predictions of Fig. lb. Setting D^i^i^z, Mp) = clearly 
has very little effect on the asymmetries, meaning that the predictions are likely to be very 
similar if NLO corrections to these contributions were included. This could be expected 
since at this cms energy, fragmentation is less the 10% of the cross section except at the 
very lowest pt values. The dotted line in Fig. 3b shows the effect of setting Rj = 1. The 
effect on the cross sections is generally less than 10% and clearly it has no effect at all on 
the asymmetries. 

Since in this calculation a jet definition is involved in NLO, the calculation of a so-called 
K-factor to estimate the size of the NLO corrections is not meaningful, as it would depend 
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on the value of Rj chosen. Thus in Fig. 3c, a comparison is made between predictions for 
the cross section using the LO and NLO matrix elements but both using the NLO structure 
functions. The aim is to examine whether the NLO matrix elements and jet definition have 
any effect on the shapes of the distributions. As stated earher, in NLO, measuring the 
kinematic variables of the photon and jet does not serve to uniquely determine the values 
of Xa and xi,. Any significant differences between the values of these variables probed in 
LO and NLO could possible show up as a difference in the shapes of the distributions. 
Figs.3a,b clearly show that there are no significant shape differences between the LO and 
NLO distributions, and therefore one can assume that estimating the values of x probed for 
a given kinematical configuration of the photon and jet by the LO formula may not lead to 
unacceptably large errors. 

In conclusion, I have examined the possibility that both the size and x-shape of the 

— * 

polarized gluon distribution of the proton, AG, may be measured at HERA-N via a mea- 
surement of the photon plus jet cross section. Control over the kinematic variables of both 
the photon and jet allows a much better determination of the a;- value probed when compared 
to inclusive prompt photon production. A comparison of the predictions obtained using dif- 
ferent polarized parton densities show that a clear distinction between both the sizes and 
shapes should be possible. Assuming that the 'fitted AG' scenario is the most plausible 
distribution, then a typical value for the asymmetry, A^^ is 10%, but given the uncertainty 
in AG the asymmetry could be as small as 2% or as large as 40%. The expected small-x 
behavior of the polarized distributions lead to predictions of negative correlations between 
the rapidities of the photon and jet. The effect is to produce very large asymmetries, even 
approaching the maximum value of ^4^,^ = 1 in certain kinematic regions, which should 
make them more easily measurable in the experiments. 
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FIGURE CAPTIONS 



[1] (a) distribution of the photon plus jet triple differential cross section 
d^a'^J IdpTpdrp'drj'^ for various polarized parton distributions for rapidities of the 
photon and jet averaged over the region —0.5 < rp ^rj"^ < 0.5. The upper dotted 
curve is for the GS A polarized distributions, (b) Longitudinal asymmetries for 
the distributions of (a), (c) and (d) similar to (a) and (b) but for both the photon 
and jet rapidities averaged over the region 0.5 < rf ,r]^ < 1.5. 

[2] (a) Distribution in the rapidity of the photon for pj. = 5 GeV and rj'^ averaged 
over the region —0.5 < r]^ < 0.5 and > pj^. (b) Asymmetries for the curves 
shown in (a), (c) and (d), similar to (a) and (b) respectively, but for 77-^ averaged 
over 0.5 <rj< 1.5. 

[3] (a) Scale dependence of the asymmetries shown in Fig. lb, for three parametriza- 
tions of the polarized densities, (b) Comparison of the asymmetries of Fig. lb 
after setting Rj — 1.0 (dotted hues) and on the same plot, after neglecting the 
contributions from fragmentation (dashed lines), (c) Comparison of the triple 
differential cross section of Fig. la with those using the LO matrix elements and 
(d) the corresponding asymmetries. 
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